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We use simple entropy arguments to estimate the possible 
size of the QGP at the AGS and the SPS. We find that the 
possibility to form a large volume of QGP at the AGS or the 
SPS is very small. The size of the QGP at RHIC and the 
LHC is also predicted. 

PACS number(s): 25.70 -q, 25.70 Pq,25.70 Gli 
One of the main goals of high energy nuclear colli- 
sions is to create a quark-gluon plasma (QGP) jl]^] of 
macroscopic size. The hope of discovering the QGP in 
heavy-ion collisions is thus to some extent connected to 
the possibility of measuring the geometric size of the re- 
gion of secondary particle production. An important tool 
in accomplishing such measurements is Hanbury Brown- 
Twiss (HBT) interferometry [|||. On the theoretical 
front, owing to the complexity of high energy heavy- ion 
collisions, several models to simulate the heavy-ion col- 
lisions process [|| @ have been devised. In those, the 
size of the QGP phase is a vital element, as it may affect 
significantly the shape of the single parton distributions. 
In this note, using HBT results from the AGS and the 
SPS, we infer the possible sizes of QGP in those energy 
regimes. The sizes of QGP at RHIC and LHC are also 
predicted. 

The idea to predict the size of QGP is simple and is 
based on the picture given by Lee some years ago [ pT"| . 
From the second law of thermodynamics we have equa- 
tion 
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Here Sqgp and Suad are the total entropy of the QGP 
phase and of the hadron phase, respectively. As more 
than eighty percent of final particles are pions, we will 
calculate the entropy of pions and multiply it by a factor 
of a to represent the total particles entropy. For this 
discussion we take a ~ 1.1. Eq.(@) can be re-written as 



sqgpVqgp < a.s v V v . 



(2) 



Here Vqqp and are the volume of the QGP phase and 
of the pion phase respectively, sqg p and s T represent the 
entropy density in the QGP phase and the pion phase, 
and reads 
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nf is the number of flavors. In the above we assume that 
the pion, quark and antiquark masses are zero. From Eq. 
(3), we have 
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Using Eq.(2) and Eq.(4), we get 

Vhad > \2VqgP- 



= 12 ~ 16. 



(4) 
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Gaussian functions have been used to fit the two-pion 
correlation functions at AGS and SPS energy and it has 
been found that the fitted source radius are 5 — 7fm || . 
From Eq.(§), we find that the source radius for QGP is 
2 — 3/rn. This result is important but depends strongly 
on the source radius as measured from pion interferom- 
etry. If the QGP size is truly 2 — 3 fm then we need to 
consider finite size effects on the spectrum distributions. 
For example, because of the Heisenberg uncertainty prin- 
ciple, particles in a small volume will have a widespread 
distribution in momentum space. It has been shown in 
Rcf. JlJ-[l^] that for quarks and gluons, finite size effects 
will become important when the size of QGP is less than 
6 fm. For pions, finite size effects set in when the size of 
the hadron phase is less than 9 fm. For such small QGP 
sizes at AGS and SPS energy one needs to consider finite 
size effects on the final state observables, and this fact is 
bound to affect several of the proposed signatures. Finite 
size effects on the observables have also been studied by 
Elze, Greiner and Rafclski |l6],[L7| long ago. If we consider 
the fact that pion interferometry actually measures the 
pion source at freezeout time and the freezeout tempera- 
ture (Tf) is less than the critical temperature (T c ), then 
the QGP source radius should be even smaller than the 
value given above. Using typical freezeout temperature 



Tf 



138 MeV or fitted results T f = 120 MeV m 



and assuming the phase transition temperature T c = 200 
MeV, we have 



V had > 37- 56 V QGP . 



(6) 



This result depends strongly on the freezeout tempera- 
ture. From Eq. (6), we find that the size of QGP phase 
is 1 — 2 fm. This conclusion implies that even if a QGP 
is formed at the AGS or the SPS, the size of QGP will 
be very small. Here we also need to mention that this 
QGP means the QGP at the critical time, i.e. when the 
system enters the mixed phase. Of course at the freeze- 
out time, the assumption that the pion mass is zero is 
inappropriate. We will come to this question at the end 
of the paper. It has been suggested that there is a pos- 
sibility to form some QGP droplets inside a hadron gas. 
If this picture is true, we have 
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Here Sdropiet is the entropy in each droplet and Nd rop iet 
is the total number of droplets. Assuming that the size 
of all droplets are the same, we have 



Vhad > (37 ~ 56) NdropletVdroplet- 



(8) 



If we assume for example that the droplet number is two, 
the droplet size will be 0.5-1.2 fm. 

The radius of pion sources at RHIC or the LHC are 
proportional to the pion multiplicity distribution as |19| : 
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From Eq.(9), we find that the hadron phase radius for 
Au + Au collisions at RHIC is 14-26 fm. First HBT 
results from RHIC indicate that the HBT radius is larger 
than the HBT radius at SPS energy, but still less than 
ten Fermis ^0|. This may due to the strong flow that 
has been observed at RHIC which make the HBT radius 
less than the true size of hadron phase, Rhad- From our 
numbers for the size of the hadron phase at RHIC we get 
an appreciable QGP size of 



Rqgp 



8/to. 
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However this estimation depends strongly on the size of 
hadron phase. This implies a caveat which we will discuss 
at the end of this paragraph. We have assumed that at 
SPS energy, the size of hadron phase is almost the same 
scale as that of the HBT radius, and this assumption is 
borne out by numerical simulations. Pion interferometry 
results from simulator (ARC, RQMD, VENUS) at SPS 
energy have indicated [glj that pion interferometry gives 
us the right geometry size of the pion source. A theoret- 
ical analysis in Ref. |2^] has shown that if we use a box 
source to fit the pion correlation function of Pb + Pb col- 
lisions at SPS energy, the box length is 12 fm, then QGP 
size at SPS energy will be 3-4 fm which is still very small. 
So both theory and empirical analyses indicate that even 
if a QGP is formed at the SPS, it's size should be very 
small and it is mandatory to consider finite size effects. 
But at RHIC energy the correspondence between physi- 
cal size and HBT radius may not hold anymore, owing to 
strong flow effects. Present HBT theory shows that HBT 
only measures part of the whole pion source due to flow 
which generates a strong correlation between coordinate 
and momentum. This causes the apparent source radius 
to become smaller. 

At LHC energy, the pion multiplicity becomes large 
|jl9f and we estimate the radius of the hadronic phase as 

Rhad = 18-36 fm, (11) 

and the corresponding QGP size should be 

Rqgp = 6-11 fm. (12) 



In the following we will use the Bjorken picture to 
perform an independent estimate of QGP size. In the 
Bjorken model, QGP is assumed to be produced at first 
and at time To it reaches equilibrium then the QGP 
will evolve according to hydrodynamical equation, that 
is MM 



s(t)t = s(t )t (r > r ). 



(13) 



When the temperature drops to T c at proper time r c , a 
mixed phase occurs, then following equation exists 



\ = (-^ (T > T ). 



s(Tc) 

According to Ref. M , after time 
T h = 6.16t c 



(14) 



(15) 



the whole system will be in a hadronic phase. This phase 
will expand and temperature will decrease from critical 
temperature (T c ) to freezeout temperature (X 1 /). During 
the hadron expansion stage, we have 



5 ( T /) T / = S ( T h)Th (T > T h ). 



(16) 



Here rt is the freezeout time. The freezeout time can be 
determined by the following equation [fLp3] 



(17) 



T ( T ) = /ZW/3 

T(r h ) t ■ 
For T f = 120 MeV and T{t h ) = 200 MeV, we have 

s{t } ) = S (r h )T h /T f = s{ Tc )C-f /3 {^p\?. (18) 

Th 1 (Th) 



Finally we get 
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(20) 



From Eq.(20) we get similar conclusion as above: even 
if there is QGP formed at AGS or SPS energy, its size 
should be around the size of nucleon. 

In the following we will use energy density of pions cal- 
culated from present measurement to estimate the pos- 
sible QGP energy density at the critical time. Using the 
Bjorken model and assuming that pion mass is zero, we 
have 



e ( T h) = ,ZLf/z 
e(r/) t h 



(21) 



Here Th and ry represent hadron time and freezeout time 
respectively. In the mixed phase, Eq.(21) is still valid, 
implying 
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Using Th = 6.16t c and assuming (see Eq. (17)) 

T ± = ( r ( r /) )3 
T f T(r h ) 

we have 



(22) 



(23) 



e(r c ) = 11.3 • £ ( T/ )(^1) 4 = 87e(r/). (24) 

When T/ = 120 MeV, it is calculated that the energy 
density of pions (tt + ,it~ and ir°) is 0.023 GeV/fm 3 . Ac- 
cording to Eq.(24), we find that the critical density of 
QGP is 



e(r c ) = 2.1GeV/fni 
Using equation (with rif = 2) 
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(25) 



(26) 



we find that the critical temperature T c = 191 MeV, 
which is consistent with the input value T c — 200 MeV. 
As we have said earlier the assumption of massless pions 
is inappropriate at the freezeout time, as the freezeout 
temperature has almost the same value as the pion mass. 
It has shown in Ref. |l^] that pion data from all heavy- 
ion reaction are consistent with thermal Bose-Einstein 
distribution 



1 



exp(£/T) - 1 



(27) 



with T = \2QMeV. Taking the pion mass as 138 MeV, 
and assuming the freezeout temperature Tf — \2QMeV , 
the entropy density for pions reads 



s = 3 



* [(! + /)!„(! + /)-/.!„/] = ^ 



(28) 



Assuming the pion mass as zero, we have s — 
0.297/(/m) 3 . That is the entropy density will decrease 
if the value of the mass increases. Due to Eq. (28), Eq. 
(0) changes to 



Vhad > 69V QGP . 



(29) 



Thus the QGP size will be even smaller than in our pre- 
vious estimates. 

Let us end with the following comments: 
(1) In the above we have assumed that the chemical 
potentials of quarks and pions are zero which is a good 
approximation if the number of quarks and pions are very 
large. From Fig. 1 or Eq. (|2^) , we find that if the chem- 
ical potential decreases the entropy density of pions will 



decrease too. When chemical potential // = m w , the 
maximum value of pion chemical potential, we need to 
consider Bose-Einstein condensate, this is of course be- 
yond the scope of the present work. From Fig.l we find 
that the corresponding maximum entropy density of pi- 
ons is around 0.6/(/to) 3 . Thus we have Vhad > 27Vqgp 
this is the upper limit of our estimation. 
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FIG. 1. Entropy density of pions, s, vs. chemical potential, 
fi. In the calculation we fixed the temperature T — 120 MeV. 

(2) In the above calculations, we have used a simple 
Bose-Einstein or Fermi-Dirac distribution for pions or 
quarks. In principle one needs to use the complete source 
distribution S(x, K) given in Ref. |2^,|24| to calculate the 
entropy density which will be used to estimate the pos- 
sible QGP size at SPS or AGS. To use the function of 
S(x, K) or f(x, K) to calculate the entropy is very im- 
portant especially when x and p correlation is essential. 
In Ref 1241, the authors have found that 



_S_ 



= 3.9 ± 1.8. 



(30) 



Taking N n = 100 - 150, we find that the QGP size is 
2 — 3/m, which is consistent with the estimations given 
above. We can use 



_S_ = Jdp[{l + f)ln(l + f)-fln(f)] 
N„ Jdpf 



(31) 



to calculate the specific entropy, S/N. Using Eq.(p7|) and 
choosing T f = 120 MeV, we find that f = 4.32. This 
result is consistent with Eq. (30). 

(3) In this paper, Eq.(|3|) has been used to calculate 
the entropy density of QGP phase in finite volume. This 
formula certainly should be corrected when QGP are con- 
fined in a finite volume. But this is precisely our conclu- 
sion: one will not observe a large extent QGP at AGS 
and SPS energies and one needs to take into account the 
finite size. In the following, we will calculate de Broglie 
wavelength of quarks, gluons and pions to show why we 
need to take into account finite size effects. Using Eq. 
@ and taking T f = 120 MeV, we find that the aver- 
age pion momentum is 372 MeV. When T c = 200 MeV, 
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the average momentum of gluons is 540 MeV and the av- 
erage momentum of quarks is 472 MeV. The de Broglie 
wavelength for pions, quarks and gluons are X n = 3.3/m, 
X g = 2.3/m and X q = 2.6/m respectively. It is interest- 
ing to notice that the de Broglie wavelength for pion at 
AGS or SPS energy is still smaller than the size of pion 
source. But where the de Broglie wavelength for quarks 
and gluons is almost the same size as the size of QGP 
source, thus finite size effects for those particles should 
be important. 

(4) In the paper, the total entropy of final state parti- 
cles is determined by the volume and entropy density of 
pions. But there are also other ways to calculate the fi- 
nal state entropy. For example, Bialas and Czyz ]2j| have 
suggested to use Renyi's entropy to derive the standard 
entropy S = — ^2 n p n ^ogp n . Here p n is the probability 
of occurrence of the state labeled by n and the sum runs 
over all states of the system. Different from the method 
given in the paper, this method does not explicitly de- 
pend on the final state volume. One can also calculate the 
final state entropy by the Landau assumption S = AN 
(consistent with Eq.(30)). Here N is the total final state 
particles. This calculation also has nothing to do with 
the final state volume. Thus, those methods are recom- 
mended to calculate the final state entropy and to infer 
the possible size of the initial state. 

(5) All the estimations in this work are based on the 
assumption that the initial state is QGP. If the initial 
state is a hot pion gas, a quark-dominated gas, or a gluon- 
dominated gas, the degrees of freedom change and the 
volume of the initial state will change too, for a given final 
state. For example for a hot pion gas, Vu na i > ^Vmitiai- 
For a quark gas, Vf ina i > 35V ini tai- For a gluon gas, 
V final > 26Vi n itai- So the initial volume is large in the 
case of a pion gas, as it should be. 

To conclude: it is a key question to know the size of 
QGP which may be produced at the early time of heavy- 
ion collisions. We show that if there is a QGP formed 
at the AGS and/or at the SPS, its radius should only 
be around 1 — 2 fm. Thus the finite size effects on the 
QGP distribution functions will be important. We also 
predict that the size of QGP at RHIC and LHC will be 
around 4 — 8 fm and 6 — 11/m respectively. For this size 
of QGP, the finite volume effects on QGP distribution 
can be neglected. 
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